Electronic structure of the MO oxides (M=Mg, Ca, Ti, V) in the GW approximation 
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The quasiparticle band structures of nonmagnetic monoxides, MO (M=Mg, Ca, Ti, and V), are 
calculated by the GW approximation. The band gap and the width of occupied oxygen 2p states in 
insulating MgO and CaO agree with experimental observation. In metallic TiO and VO, conduction 
bands originated from metal 3d states become narrower. Then the partial densities of transition 
metal Cg and t2g states show an enhanced dip between the two. The effects of static screening 
and dynamical correlation are discussed in detail in comparison with the results of the Hartree- 
Fock approximation and the static Coulomb hole plus screened exchange approximation. The d-d 
Coulomb interaction is shown to be very much reduced by on-site and off-site d-electron screening in 
TiO and VO. The dielectric function and the energy loss spectrum are also presented and discussed 
in detail. 

PACS numbers: 71.10.-w, 71.15.-m, 71.20.Be 



I. INTRODUCTION 

The first-principle electronic structure calculation, 
based on the local-density approximation (LDA) within 
the density- functional theory (DFT) ,EJ has had great suc- 
cess in theoretical investigations of the ground state prop- 
erty of condensed matters. However, in many LDA cal- 
culations, the band gap is unsatisfactory underestimated 
in insulators and semiconductors. In the electronic struc- 
ture of 3d transition metals calculated by the LDA, the 
occupied 3c? band width is too broad, the exchange split- 
ting is overestimated and the satellite structure in the x- 
ray photoemission spectroscopy (XPS) does not appear. 
These quantities are associated with quasiparticle prop- 
erties. The transition metal monoxides TO (T=Mn, Fe, 
Co, and Ni) are wide-gap antiferromagnetic insulators 
in experiments, but FeO and CoO appear as metals, and 
MnO and NiO as small-gap semiconductors in the LDA.cl 
The magnetic moments of TO are underestimated by the 
LDA. 

The electronic structure in strongly correlated elec- 
tron systems has been highly challenging field. The self- 
Coulomb and self-exchange interactions almost cancel 
with each other inj-piany cases of the LDA but still not 
completely in TO a The self- interaction corrected LDA 
(SIG-iLDA) has been applied to transition metal monox- 
idesD^El It gives satisfactory results for the band gap and 
magnetic moments in the transition metal monoxides 
TO (TrgjMn, Fe, Co, and Ni) except VO. The LDA-HU 
methodQCla works reasonably well for the Mott-Hubbard 
insulators or rare earth metal compounds where the id 
or 4/ bands are partially filled. The LDA-I-U method is a 
static limit of the GW approximation (GWA).la However, 
application of the LDA-I-U method to early transition 
metal monoxides TiO and VO would lead to antiferm- 
magnetic insulators rather than paramagnetic metalsE 

The GWA is the first term approximation of the many- 
body perturbation series for the Green function, and is a 
simple but excellent method to leakulate the self-energy 
with the dynamical correlation.clll2l On the other hand. 



the HF approximation overestimates the band gap in in- 
sulators and semiconductors because of the lack of the 
screening terms and gives zero density of states at the 
Fermi level in the electron gas. In the GWA, there is no 
self-interaction because the exchange term of the GWA 
is of the Hartree-Fock type. The dynamical correlation 
effect in the GWA is treated by the random-phase ap- 
proximation (RPA). In the GWA, the dynamical effects 
give rise to important features that are not accounted 
for by the HF approximation, for example, damping of 
quasiparticles and characteristic structure in the spectral 
function. The screening effects recover the finite density 
of states at the Fermi level. The band gap is corrected 
and reduced in the GWA by the dynamical screening 
effects in comparison with the one of the HF approxima- 
tion. The dynamical correlation effect is, important for 
the band width in the transition metals U We could not 
expect improvement in the GWA for the exchange split- 
ting and satellite structure because it needs higher order 
diagrams (e.g., vertex correctionsl-for electron-electron 
and hole-hole scattering proccsses.Ej 

The GWA has been applied to several real systems, 
for f^ypmplp simple metalsji3 semiconductors and insu- 
latorsJlJIlj' !^ transition metalsjl3 and transition metal 
monoxides£3 In simple metals and semiconductors, the 
GWA can be formulated with plane wave basis set based 
on the pseudopotential method and the plasmon pole ap- 
proximation lias been often used for calculating dielec- 
tric function.tj In contrast, transition metals and their 
compounds have strong atomic potentials and the GWA 
based on the pseudopotential method cannot be applied 
because of restriction of the number of plane waves. It is 
essentially important to include core electrons in many 
cases and, moreover, the plasmon pole approximation is 
not applicable to the dielectric function for transition 
metals and their compounds because there is no isolated 
plasmon pole. Therefore, all-electron calculation and no 
plasmon-pole approximation in the dielectric function are 
required for transition metals and their compounds. 

The GW calculation of Ni, based on the linearized 
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augmented-plane-wave (LAPW) basis and no approxima- 
tion in the dielectric function, gives results in good agree- 
ment with experiments.E-O The modeling of the screened 
interaction in the GWA is able-to explain satisfactorily 
photocmission spectra in MnO|13 NiOj23 and V02-£il In 
the pioneering work of the GW calculation of NiO by 
Aryasetiawan, the band gap and magnetic moment agree 
quite well with experimental values.Ej In the analysis of 
photocmission spectroscopy by the cluster configuration- 
interaction theory, NiO is assigned as a charge-transfer 
type insulator ,E3l but in the LDA it is classified as a Mott- 
Hubbard type insulator because the transition metal 3d 
bands appear at higher energy than oxygen 2p bands .□ 
The relative position of the oxygen 2p states of the GWA 
is also improved in comparison with that of the LDA. 
However, the top of the valence bands in the GWA is 
mainly Ni id states and NiO could not be a charge- 
transfer type insulator. 

In most GW calculations, full self-consistency is not 
carried out. The self-consistency in the GW calculations 
can maintain the conseipffition laws of particle number, 
energy, and momentumfij However such treatment can- 
not gttarantee better agreement with experimental re- 
sultsB 

Because it is derived from many-body perturbation 
theory, the GWA is applicable to wide variety of classes 
of systems, and should be tested in various systems. In 
this paper, the GW methodj-based on the linear muffin- 
tin orbital (LMTO) methodO is applied to several non- 
magnetic oxides, insulating MgO, CaO and metallic TiO, 
VO. Only the electronic structure in MgO was calculated 
by the GWA beforeJl3 The insulators MgO and CaO are 
typical simple oxides, whose insulating band gap is un- 
derestimated by the LDA. TiO and VO are nonmagnetic 
metallic oxides. In the LDA their characteristics are very 
similar. We will discuss the electronic band structure and 
its systematic variation for MgO, CaO, TiO, and VO, 
concerning with the iJ-k relation, band gap, band width, 
the density of states, plasmon frequencies, structures in 
the energy loss spectrum, and other properties. We will 
also discuss the correlation effects, including the strength 
of the Coulomb interaction with d-d screening effects. 

The paper is organized as follows. A simple review on 
the theoretical framework will be given in Sec. |l| The 
results for these systems are presented and discussed in 
Sec. HI. The dielectric function and the energy loss spec- 
trum are also discussed. Finally, in Sec. IV we present 
our summary. 



II. THEORY 
A. GW approximation 

When effects of dynamical correlation are included, the 
self-energy is expressed formally as a series expansion of 
dynamically screened interaction. The GWA replaces the 
self-energy by the lowest order term of the expansion as 



W{1,2)= / d(3)e-^(l,3)i;(3,2). 



E(l,2) = iG(l,2)W(l,2). G is the one particle Green 
function and the dynamically screened interaction W is 
defined by 

(1) 

= «(1,2) + y d(34)i;(l,3)x°(3, 4)1^(4,2), (2) 
= «(1,2) + T^^(1,2), (3) 

where e is the dynamical dielectric function, v is the 
bare Coulomb potential and is the irreducible po- 
larization function x°( 1,2) = -iG(l, 2)0(2, 1). Here 
we use an abbreviated notation (1) = (ri,(Ti,ii) and 
w(l,2) — v{vi,V2)5{ti — t2). Equation (||) can be treated 
by the RPA. 

We adopt the LDA Hamiltonian to be the unperturbed 

one H'^ , 



H° = T + V 



H 



LDA- 



(4) 



Here T is the kinetic energy, is the Hartree poten- 
tial, and V|^DA is the exchange-correlation potential in 
the LDA. We presume the wave functions {'0kn(r)} of 
the LDA to be a reasonably good starting wave functions, 
though this assumption should be carefully investigated 
in detail. Then the self-energy can be written by three 
terms as 



AS = 1]"= + - K 



LDA, 



(5) 



where E^(= iGv) is the exchange part (the Fock term) 
and E'^(= iGW^) is the dynamical correlation part. The 
quasiparticle energy is given as 



Ek7i — Ekn + ^knAEkn(ekn) 



(6) 



where e^n is the LDA eigenvalue. The self-energy AEkn 
and the renormalization factor Z\^„, the weight of the 
spectral function, are given by 

AEk„(ek„) - (^^k„|E" + E^(ek„) - l^^AlV^kn), (7) 



9AEk„(c^) 



duj 



(8) 



The renormalization factor Zkn is a measure of the oc- 
cupation number and should equal to the discontinuity 
of occupation number at the Fermi energy. Therefore 
it should satisfy the condition Z^n ^1- In the present 
work we perform one iteration calculation without self- 
consistency. 

In order to know the difference of the dynamical and 
static screening effects, we also apply the static Coulomb 
hole plus screened exchange j^COHSEX) approximation 
for the self-energy by Hedin.l3 The static COHSEX ap- 
proximation can be obtained by neglecting energy depen- 
dence of the self-energy. Then the self-energy is expressed 
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TABLE L The band gap Eg, the width of valence band Wv, and the position of the O s band Eq, 
CaO. The unit is (eV). 



in insulating MgO and 





LDA 


HF 


MgO 

COHSEX 


GWA 


expt. 


LDA 


HF 


CaO 
COHSEX 


GWA 


expt. 


Eg 


5.2 


18.4," 16.5," 17.6" 


9.6 


8.2 


7.8" 


3.65 


15.9" 


7.7 


6.64 


7.1" 




5.0 


10.4," 5.8," 7.64" 


5.6 


5.0 


5.0-6.0'* 


2.8 


3.43" 


3.4 


2.9 






~ 16 


~ 40," - 24" 


~ 19 


~ 17 


18-21" 


~ 15 




~ 17 


~ 16 





"Reference 2t . 

"Reference 27 . I . 

" (Uncorrelated) HF results in Ref. GSl Correlated HF results are 
Eg=S.21 eV in MgO and _Bg=7.74 eV in CaO. 
"^Reference 2a. 
"Reference 30. 



as 

occ 

SsEx(r,r') = -^^k„(r)Vk„(r')W^(r,r';^^ = 0), (9) 

kn 

ScOH(r,r') = i5(r-r') [t^(r, r'; = 0) - ^;(r, r')] , 

(10) 

where Ssex is the screened-exchange term and Scoh is 
the Coulomb-hole term which becomes the local potential 
in the static approximation. The self-energy matrix ele- 
ment of W{r, r' : w = 0) in the d orbitals is the Coulomb 
interaction including the d-d screening, which will be dis- 
cussed in the following section. 



B. LMTO minimal basis set and choice of LMTO 
parameters 

Because the plane wave basis becomes very costly 
for sysijems containing 3d or 4/ electrons, the LMTO 
methodc3 is more appropriate. We use the LMTO ba- 
sis set XRLj/(r) within the atomic sphere approximation 
(ASA) for the LDA calculation. Here L is angular mo- 
mentum L = (/, m). The LMTO can be expanded by 
the mufhn-tin orbital 0RLi/(r) and it's energy derivative 
^RL!/(r)- The product basis formahsmO is used and 
proved to be very efficient. Computational details and 
numerical technique are shown in Appendixes ^ and 

The lattice constants of MgO, CaO, TiO, and VO axe 
a = 4.2122, 4.8105, 4.1766, and 4.062 A, respectively.El 
Empty spheres (ES's) are often necessarily introduced 
in the ASA formalism except in closed packing lattices. 
The calculation of insulating MgO and CaO without ES 
gives a larger GW band gap by about 1 eV than that 
of calculation with one ES. In the metallic TiO and VO, 
proper states arc not obtained by the calculation without 
ES's in the LDA. In the case of two or more ES's in NaCl 
type structure, the radii of ES become too small. In the 
calculation with two ES's the band gap is almost the 
same as in the calculation with one ES. Therefore, the 
one type of ES is enough in the NaCl type structure. In 



the present paper, we use the following sets of radii (in 
atomic unit) of atomic spheres; (Mg, O, ES)=(2.21, 2.21, 
1.62), (Ca, O, ES)=(3.03, 2.19, 1.49), (Ti, O, ES)=(2.71, 
1.87, 1.13), and (V, O, ES)=(2.58, 1.91, 1.12). In metallic 
systems we use 512 (8x8x8) mesh points to have enough 
convergence. On the other hand, in insulating systems 
smaller number of mesh points 64 (4 x 4 x 4) is good 
enough to achieve a convergence because of the existence 
of the band gap. 

We must treat the wide energy range in the calculation 
of the dielectric function and, therefore, use a sufficient 
number of LMTO's. The energy-dependent polarization 
function x^{^) h^^s a long tail in higher energy range 
and, therefore, metal / and oxygen d or f orbitals are 
included for the GWA in the present work, though these 
orbitals are omitted usually in the LDA. The existence 
of ES is important but higher order orbitals in ES are 
not essential. The set of the maximum orbitals I of the 
LMTO basis in M, O and ES atomic spheres are chosen 
to be (//s) for MgO, (ffd) for CaO, (fds) for TiO, 
and {fds) for VO. In CaO, since the band gap becomes 
narrower, we include the basis of ES with higher angular 
momentum. There are no essential differences between 
the different choices but the above choice is the best for 
the band gap within ±0.1 eV. 



III. RESULTS AND DISCUSSIONS 
A. Energy band structure 

1. Insulator : MgO and CaO 

MgO and CaO are ionic crystals and are band insu- 
lators. In MgO, the band gap separates the unoccu- 
pied oxygen 2p valence bands and the unoccupied mag- 
nesium 3s-3p conduction bands. In CaO the band gap 
locates between the unoccupied 4s-4p conduction bands 
and the unoccupied calcium ?>d conduction bands. The 
band structures of MgO and CaO, calculated both in 
the LDA and the GWA, are shown in Fig. Q along high 
symmetry lines. The band gap, the width of the upper 
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FIG. 1: The energy bands of MgO (left) and CaO (right), calculated by the LDA (solid lines) and the GWA (dashed lines) along 
high symmetry lines. The closed circles are the calculated points in the GWA. The high symmetry points are L=(l/2, 1/2, 1/2), 
r=(0, 0, 0), X={1, 0, 0), W={1, 1/2, 0), and A"=(3/4, 3/4, 0). 



valence bands and the position of the low lying oxygen 
2s band are summarized in Table |, in comparison with 
those by the LDA and the HF approximation. Our re- 
sults of MgO are in good-agreement with those of the 
previous GW calculation.ll3 

The band gap in the GWA is found to be 8.2 eV in 
MgO and 6.64 eV in CaO, which are in good agreement 
with the experimental values 7.8 and 7.1 eV,E3 respec- 
tively. The band gap in MgO is the direct one at F point 
both in the LDA and the GWA as in experiments. The 
band gap in CaO is the indirect one between F and X 
points in the GWA as in the LDA, although the energy 
difference between the F-F gap and T-X gap becomes 
smaller in the GWA. The conduction band minimum at 
the X point originates from metallic d states and the 
lowering the conduction band minimum at the X point 
is due to deepening of the potential of the metal ion. On 
the other hand, the direct gap of CaO at the F point was 
proposed by the HF approximation with correlaticwi cal- 
culated by the second order perturbation theory.c3 We 
are not aware of any angle-resolved photoemission or in- 
verse photoemission experiments on CaO, and the nature 
of the fundamental gap of CaO cannot be concluded at 
present. 

The GWA valence band width is 5.0 eV in MgO and 
2.9 eV in CaO, which are in good agreement with ex- 
periments and almost the same as those of the LDA. 
The LDA valence band width is generally narrower in 



general than the GWA and the observation, because the 
wave functions of valence states are more concentrated 
on atomic sites by local approximation in the LDA, and 
the effect e£|the exchange and correlation effects are much 
enhanced.E3 However, because the valence bands are oxy- 
gen 2p states, such exchange and correlation effects do 
not affect much in MgO and CaO and band width re- 
duction is not much. In the-calculation with the exact- 
exchange (EXX) potentialfJ the valence band width of 
MgO and CaO is narrower than the one in the LDA but 
this is due to the different reason. In the EXX calcula- 
tion, the exchange energy is approximated by the spher- 
ically symmetrized Fock term and the correlation energy 
is replaced by that of the LDA. The EXX potential is 
larger in comparison with the LDA exchange potential 
and gives narrower valence bands of MgO and CaO. In 
the HF approximation, which includes the non-local in- 
teraction, the valepce band width is broadened, for in- 
stance, in MgOcj'EZl But the band gap and band width 
are overestimated in the HF approximation due to the 
lack of screening effects. In the GWA the valence band 
width is reduced in comparison with the HF method. 

The position of oxygen 2s band is pushed down in the 
GWA so it is close to the experimental value. The LDA 
makes the core or semicore states locate into shallow en- 
ergy regions, compared with experiments, because the 
LDA potential contains self-interaction. 

The self-energy and the correction ZAY, to the LDA 
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TABLE II: The band width of O 2p band W2p, the position of O s band Eq^ , and the band width of transition metal 3d band 
Wsd, Wt2g, and W^^ in metaUic TiO and VO. The unit is (eV). 







TiO 






VO 






LDA 


COHSEX 


GWA 


LDA 


COHSEX 


GWA 




4.6 


4.7 


4.1 


5.2 


5.3 


4.6 




~ 22 


~ 25 


~ 23 


~ 22.5 


~ 25 


~ 23 


W3d 


8.4 


9.1 


6.8 


8.1 


7.6 


6.4 




6.5 


6.1 


4.3 


6.4 


5.4 


3.8 


We, 


4.1 


4.1 


3.2 


4.5 


4.4 


3.3 



energy of MgO are shown in the left side of Fig. g. The 
gap at the top of the valence bands (the energy zeroth) 
is seen in the self-energy and as a result the energy cor- 
rection ZAE has the discontinuity. In conduction bands, 
fluctuations of and V|^da caused by different or- 
bital characters, but the correction terms ZAS converge 
to almost constant value. The correction terms of the 
highest valence band Zhvb and those of the lowest con- 
duction band .^lcb are as follows: in MgO, Zhvb = 0.75- 
0.8 and Zlcb = 0.83-0.87; in CaO, Zhvb = 0.73-0.75 
and Zlcb = 0.80-0.85. The smaller Z of CaO shows 
that the electron-electron correlation in CaO is stronger 
than in MgO. 

To see the dynamical screening effects in MgO and 
CaO, we also calculated the energy band structure by 
using the static COHSEX approximation and shown in 
Table |. In the static COHSEX approximation, the band 
gap is much smaller than the HF results but is still over- 
estimated to be 9.6 eV in MgO and 7.7 eV in CaO. The 
reduction of the band gap from the HF result to the static 
COHSEX result and that from the static COHSEX result 
to the GW result are attributed to the static screening 
effect and the dynamical screening effect, respectively. 
Therefore, in MgO and CaO, the static screening is the 
major contribution rather than the dynamical effects. 
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FIG. 2: Energy dependence of the self-energy of MgO (left) 
and TiO (right). The closed diamond for VldA) the open 
triangles for E^, the open squares for E"^, and the closed circles 
for ZAE. 



The valence band width is 5.6 eV in MgO, and 3.4 eV 
in CaO in the static COHSEX approximation. In the 
GWA the valence band width is slightly reduced in com- 
parison with the static COHSEX approximation. The 
reduction from that of the HF approximation to that 
of the static COHSEX approximation is much larger in 
MgO, because of the large difference of electronegativity 
and resultant importance of the static screening. The 
comparison of the band width by the GWA with that by 
the static COHSEX approximation clearly shows that the 
dynamical effects are appreciable in the electron-electron 
correlation. 



2. Metal : TiO and VO 

TiO and VO are metals, showing typical properties of 
covalent crystals such as high melting temperature and 
extreme hardness. The band structures of TiO and VO 
are shown in Fig. ||, calculated both in the LDA and 
the GWA and summarized in Table ||. In the LDA the 
oxygen 2s states are located around about 22 eV below 
Ep in TiO and VO. These oxygen 2s states are separated 
by about 11 eV from the bonding states consisting of 
metal id and oxygen 2p orbitals in both TiO and VO, 
mainly oxygen 2p, at 7 ^ 10 eV below Ep- A real gap 
of about 2 eV in TiO and that of about 1.5 eV in VO 
can be seen between this bonding states (mainly oxygen 
2p) and the antibonding states (mainly transition metal 
3d) at around Ep- The occupied conduction states are 
heavily dominated by transition metal 2>d orbitals. 

The band width of metal Sd states becomes narrower in 
the GWA than in the LDA. The occupied id band-width 
from the UPS and XPS is about 3 eV in TiO1.03.Ea The 
width of the unoccupiedrSd band in the GWA agrees with 
that of the BIS spectra.Ea 

In TiO and VO the width of oxygen 2p valence band 
in the GWA becomes narrower in comparison with that 
of the LDA, which is not the case in MgO and CaO. The 
hybridization is very small in MgO and CaO. In TiO and 
VO, the hybridization between the transition metal Zd 
and oxygen 2p increases due to deepening metal 3d states 
with the increase of the d electron. Once one adopts the 
GWA, the wave functions of valence bands are localized 
more due to the screening effects and the width of the 
band around —10 eV becomes narrower. Then the con- 
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FIG. 3: The energy bands of TiO (left) and VO (right), calculated by the LDA (solid lines) and by the GWA (dashed hnes) 
along high symmetry lines. The closed circles are the calculated points in the GWA. The highest energy band in the figure is 
mainly from metal s states. For the symmetry points, see the caption of Fig. hi. 



duction bands are lifted to higher energy side, relative to 
the lower energy oxygen 2p (bonding) band in TiO and 
VO. The effects of the reduction of the conduction (tran- 
sition metal 3d) band width also causes the increasing 
the separation between eg and t2g states. The position 
of oxygen 2s band is also pushed down by the GW correc- 
tion, and this is caused by the self-interaction correction 
as in MgO and CaO. 

The self-energy and the correction ZAY, to the LDA 
energy of TiO are shown in the right side of Fig. ^. In 
contrast to insulating system MgO three parts of the self- 
energy and the energy correction ZAS are continuous at 
the Fermi energy Ep, which are the typical properties of 
the Fermi liquid. The renormalization factor of transi- 
tion metal 3d states is Z = 0.55-0.65 in TiO and Z — 
0.51-0.63 in VO. Z of oxygen 2p states is about 0.7-0.8 
in both systems. Those results of the renormalization 
factor suggest that the interaction between 3d electrons 
is strong, and the correlation in VO is larger than in TiO. 

The d-d Coulomb interaction {(fidipdlW {ui = 0)\(f)d4'd) 
is calculated to be about 1.5 eV both in TiO and VO. 
The bare Coulomb interaction {(t>d4>d\'v\4>d4'd) is actually 
15.3 eV in TiO and 17.8 cV in VO. Then the correla- 
tion term {cjidM^i^ = 0)\Md) is -13.7 eV in TiO 
and —16.3 eV in VO. Therefore, the correlation effects 
and the static screening (at the Fermi energy) are quite 
important in these metallic oxides. We should mention 
that the d-d Coulomb interaction {(j)d4>d\W {u; = 0)\4>d4'd) 



is different from the Hubbard U evaluated from the LDA, 
which includes only screening by on-site d electrons. The 
term {(j)d(j)d\W {uj — 0)\(j)d4'd) includes the screening ef- 




Energy (eV) Energy (eV) 



FIG. 4: Total and projected DOS in a transition metal atom 
in TiO and VO by the LDA (bottom) and the GWA (top). 
The solid lines for the total DOS in a unit cell, the dotted 
lines for the partial DOS of transition metal Sd-Cg states and 
the chain lines for the partial DOS of transition metal 3d-t2g 
states. The contribution by the oxygen atom is small in this 
energy range. 
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FIG. 5: From left to right, spectra of MgO, CaO, TiO, and VO are shown. Upper panel and lower panel are the dielectric 
function e{uj) and the electron energy-loss spectra (EELS) e~^{ui) for q — 27r/a(0.25, 0, 0). 



fects by both on-site and ofF-site d electrons.E3 

The density of states (DOS) of transition metal 3d 
orbitals are shown in Fig. ^. In the calculation of the 
density of states in the present GWA, the eigenvalues 
were replaced by those of the GWA but the eigenfunc- 
tions were the same as those of the LDA, and the lifetime 
(from imaginary part of the self-energy) of quasiparticle 
was not considered. The transition metal eg states are 
pushed up in the GWA and the eg-t2g separation becomes 
much clearer than in LDA results. 

To see the dynamical correlation effects in TiO and VO 
we also show the results by the static COHSEX approx- 
imation in Tabic The band width of the static COH- 
SEX approximation docs not become narrower apprecia- 
bly in comparison with the LDA results. The reduction 
by the GWA is much appreciable. We could conclude 
that the band width reduction is mainly due to the dy- 
namical correlations. The importance of the dynamical 
correlation was also mentioned in thje 3c? band width of 
3d ferromagnetic transition metals.tll 

Apphcation of the LDA-I-U method to TiO and VO 
lead to antiferromagnetic insulators rather than param- 
agnetic metals, when one uses the U vaJue obtained by 
the LDA and the photoemission spectrafl However, as al- 
ready shown, the d-d Coulomb interaction {(j>ii(pii\W'^{uj = 
O)|0d0d) is much smaller by the correlation effects. 



B. Energy loss spectrum and plasmon frequency 

The dielectric function is seen in the procedure of the 
GW calculation. We show the dielectric functions and 
the electron energy-loss spectra (EELS) in Fig. || for q = 
27r/a(0.25, 0, 0), using 512 k points in the Brillouin zone. 

The dielectric function consists of a part of free elec=. 
tron contribution and a part of interband transition.EZI 
The plasmon frequency Wpi corresponds to the zeroth of 
the free electron part of the dielectric function and the 
peak of the EELS. When the interband transition ex- 
ists at lower energy side of the plasmon frequency, the 
frequency of vanishing £(q, w) shifts towards higher fre- 
quency side. This is the shift of the plasma edge in the 
dielectric functions of metals or the shift of the plasmon 
peak in the EELS. The peaks of the imaginary part of 
the dielectric function e{uj) appear essentially at the en- 
ergies of interband excitations. On the other hand, the 
peaks of the EELS appear at positions where the dielec- 
tric function vanishes or becomes very small. Therefore, 
the EELS have complex structure originated from the sit- 
uation of the coexistence, in a narrow energy region, of 
the plasma edge and the interband transition. This is the 
case in transition metal monoxides and it may be very 
interesting to see the difference or the similarity of the 
dielectric functions and the EELS of MgO, CaO, TiO, 
and VO. 

The interband transition from O 2p to metal 3d ap- 
pears at about 10 eV in CaO, TiO, and VO but that of 
MgO locates at much higher energy. In MgO, the tran- 
sition from O 2p to metal 3s-3p creates a broad peak of 
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e{u) at about 12 eV. On the other hand, the transition 
from O 2p to metal 4s-4p makes a structure of e{u!) at 
about 15 « 20 eV in CaO, at about 20 « 27 eV in TiO 
and VO. The sharp peaks of e(w) at 1.5 eV in TiO and 
at 1.1 eV in VO are created by metal 3d-3d transition. 
The dielectric function of metals should reduce generally 
to zero near lo = and then show the Drude peak at 
w = 0. In metallic compounds TiO and VO, low energy 
excitation peaks of e{ui) exist down to 0.2 « 0.3 eV. Then 
the nonzero value of lme{uj = 0) in Fig. || is the artifact 
caused by broadening these low frequency peaks of the 
irreducible polarization function by Gaussian with a 
width of a w 0.7 eV. 

In the EELS, the sharp peak of metallic TiO and VO 
at 5 eV is originated from metal 3d-3c? transition. The 
transition from O 2p to metal 3d generates the peak at 
14 eV in CaO and at 18 eV in TiO and VO. The shoul- 
der at around 17 eV in MgO is caused by the transition 
from O 2p to metal 3s-3p because the imaginary part of 
dielectric function has a relatively large nonzero value at 
this energy due to broad peak of interband transition. 
The plasmon peak shifts from 21.1 to 25.0 eV in MgO, 
from 17.2 to 21.8 eV in CaO, from 24.6 to 27.8 eV in 
TiO, and from 27.2 to 29.9 eV in VO due to interband 
transitions. In this estimation of the plasmon frequency, 
we adopted the number of electrons to be 6 for MgO and 
CaO (metal s and O 2p electrons), 8 for TiO (Ti 3d,4s 
and O 2p), and 9 for VO (V 3d,4s and O 2p). 



IV. SUMMARY 

We applied the GWA to the nonmagnetic oxides MgO, 
CaO, TiO, and VO. Band gap and occupied oxygen 
2p band width in insulating MgO and CaO are in 
good agreement with experimental results. In both sys- 
tems oxygen 2s level is pushed down, because the self- 
interaction is corrected in the GWA. We compare the 
results by the GWA with those by the static COHSEX 
approximation and discussed the effects of the dynamical 
correlation. In metallic TiO and VO, the band width of 
transition metal Cg and <2g states become narrower due 
to the correlation effects. 

We also present values of the d-d Coulomb interaction 
in TiO and VO and showed the crucial role of the dynam- 
ical correlation in these systems. The GW method may 
be one of the most appropriate method which can pro- 
vide the information of the electron-electron interaction 
with screening effects. 

The actual wave functions should differ from those in 
the LDA when the effects of dynamical correlation are 
included. The choice of the starting wave functions or 
equivalently the choice of the unperturbed Hamiltonian 
in the GWA may be very important from this point. This 
may be very essential for the characteristics in the mag- 
netic transition metal oxides but still open for the future 
study. 
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APPENDIX A: COMPUTATIONAL DETAILS 

Since the wave function ip is written in terms of the 
LMTO basis functions x, the Hilbert space of the self- 
energy is spanned by the product basis as 

{S} = {^A^} = {xx} = {00} + im + {00} • (Al) 

In the present calculations, because of a proper choice of 
the energy parameters the absolute values of the mixing 
coefficients of to are less than 0.1 for all terms and 
the norm of 0, (00), is 0.1 ~ 0.3 even in the largest case. 
Then the terms including can be dropped out safely and 
the self-energy E can be expanded by 0;*„,^, (r)0/m^(r) 
= 0/v'(|r|)0i,.(|r|)y;,*„,(f)y,„(f). The product of two 
spherical harmonics is reduced to the linear combination 
of one spherical harmonics. Consequently we can con- 
struct the new basis set 

-^R.Z'!v'!y;Z"m"(''R.) 

= ^e*-T0z.,.(|rRT|)0;.(|rRT|)i1",n"(f5S), (A2) 

T 

where L" is depend on I' and I, that is \l — 1'\ < I" {I', I) < 
I + I' and \m"{l\l)\ < and trt = r - R - T. 

The orthonormal basis set B can be produced by the 
linear combination of the nonorthogonal basis set B as 
Bf = J2jBfD^i, where j = (R/VZz/; Z"m"), i is a new 
sufhx in the orthonormal basis set and D is the trans- 
formation matrix from the nonorthonormal basis to the 
orthonormal basis. Here we use the convenient notation 

\Bf) ^ - |^k)^^ Kk> ^ \Bf). (A3) 

The matrix is obtained by orthogonalizing the overlap 
matrix Ofj = («kbk)- 

We can expand the self-energy in terms of the new 
orthonormal basis set |ik). To calculate the self-energy 
we first calculate the polarization function with the 
nonorthonormal basis. 

The Coulomb potential is calculated by using the 
bare structure constant of the LMTO method.E3 In this 
method, two parts of the Coulomb matrix can be calcu- 
lated separately, on-site and off-site parts, by using the 
LMTO bare structure constant S'°(q). 

Next we need to transform from nonorthonormal basis 
|zk) to orthonormal basis |zk), that is (i|x'^|j) — > (*lx"b') 
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and — > {i\v\j). Within the orthonormal basis we 
can calculate the response function X = (1 ~ X°'f^)~^X° 
and the screened Coulomb potential — vxv. Finally 
the self-energy can be obtained by the inverse transfor- 
mation from the orthogonal basis to the nonorthogonal 
basis, and that from the product basis to the Bloch basis. 

The number of basis set can be reduced by excluding 
terms of higher I" . We have studied the dependence on 
a choice of the maximum l"'s for set of I" in (Mg, O, ES) 
atomic spheres being (3,3,0), (4,4,0) and (6,6,0). The 
configuration (6,6,0) corresponds to the full calculation 
for orbital / for metals, / for oxygen and s for ES. The 
total number of product bases of (3,3,0), (4,4,0), and 
(6,6,0) are 120, 174, and 222, respectively. The proper 
choice may be a set of the maximum Z"'s which produces 
the total number of the product basis nearly equals to a 
half of that in the maximum choice. For example, the 
product basis of (3,3,0) may be enough for the set of 
the ffs Bloch orbitals. In present work the full product 
basis is selected, but this way of product basis reduction 
is very effective and powerful for large systems. 

APPENDIX B: NUMERICAL TECHNIQUE 

In the calculation of the self-energy, the integration is 
performed over the whole Brillouin zone. The Coulomb 



matrix w(q) has a singularity at q = as F{q) = l/|qp. 
The integration of f (q) over the Brillouin zone docs not 
diverge but special cares are needed not only for the 
q = term but for small finite q. In the present calcu- 
lation, the integration over the Brillouin zone is replaced 
by the summation over discrete points. We use a set of 
discrete points distributed densely near the q = point, 
not on uniform k mesh, in the Brillouin zone. For a 
choice of the discret^JDoints near q = 0, we use the off- 
set F-point method.Ej The integration of F{q) over the 
Brillouin zone can be performed analytically and the off- 
setted points Q's are chosen near q = so as to satisfy 
a relation 

/ F(q)dq = ^F(Q) + ^F(k), (Bl) 

■^B^ Q k^a 

where k's are the discrete mesh points of the Brillouin 
zone. The choice of the offsetted points Q are taken into 
consideration so as to keep the symmetry of the system. 

We tested the offset method in the exchange energy 
of the electron gas system and the accuracy has been 
confirmed even in case of small number of k-mesh points. 
The careful treatment of the Coulomb matrix at or near 
q = is very crucial near the band gap or the Fermi 
level. 
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